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SOME NUMERICAL RESULTS ON THE BEHAVIOR OF 
ZEROS OF THE HERMITE-PADE POLYNOMIALS 


N. R. IKONOMOV, R. K. KOVACHEVA, AND S. P. SUETIN 


Abstract. We introduce and analyze some numerical results obtained 
by the authors experimentally. These experiments are related to the 
well known problem about the distribution of the zeros of Hermite- 
Pade polynomials for a collection of three functions [fo = 1, fi, T 2 ]- The 
numerical results refer to two cases: a pair of functions /i,/2 forms 
an Angelesco system (see (3)) and a pair of functions /1 = f,f 2 = / 2 
forms a (generalized) Nikishin system (see (4)). The authors hope that 
the obtained numerical results will set up a new conjectures about the 
limiting distribution of the zeros of Hermite-Pade polynomials. 
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1. Introduction 


1.1. In the present work we introduce and analyze some numerical results 
obtained by experimenting. These experiments are connected with the well 
known problem about the distribution of the zeros of Hermite-Pade poly¬ 
nomials for a collection of three functions [fo = 1, / 1 , jg], defined and holo- 
morphic at the infinity point z = 00 , / 1 , /2 £ J^(oo). Our numerical exper¬ 
iments are restricted only to the Hermite-Pade polynomials of first kind. It 
is well known that Hermite-Pade polynomials 1 of first and second kind are 
closely related to each other, see [46, §2, equation (2.1.9)], [13], [14], [15]. 
The numerical results obtained here may permit interpretations about the 
Hermite-Pade polynomials of the second kind. 
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Sometimes these polynomials are called “multiple orthogonal polynomials”, see [14], 
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Let P n := C n [z\, n E N, be the class of all polynomials with complex 
coefficients of degree ^ n. For an arbitrary n E N, define the Hermite- 


Pade polynomials of first kind Q 


n,j 


n j — Oj I; 2, Q n , i) Qn, 2 ^ 0, by the 


relation (see [46], also [24], [60], [28], [2], [8], [3], [14]) 
( 1 ) ( Qn,0 ■ 1 + Qn,lfl + Qn,2f2){z) = O ^ y2n+2 


oo. 


In the present work, we will use the terminology and notation of the work 
by J. Nuttall [46] (see also H. Stahl [60]). According to their paper, the 
classical Pade polynomials P n ,o,Pn,i of the function / E J^f(oo) are, in 
fact, the Hermite-Pade polynomials for the collection of two functions [/o = 
1, fi = /], such that Pn, o, P n ,i G Pn, P nj i ^ 0 and: 

(2) ( P n, 0 • 1 + Pn,lf){z) = O , Z -Aoo. 


We present the results from our numerical experiments about two cases. 
In the first case the functions f\ and f 2 have the following form: 


( 3 ) 


fi{z) 


z 

\j{z- a\)(z — 6i) ’ 


h{z) 


z 

Viz- a 2 )(z - b 2 y 


where a\, &i,a 2 ,&2 E C, ai / b\, a 2 ^ b 2l and {ai, b\} n {a 2 , b 2 } = 0- 
Therefore, the pair of functions fi, f 2 forms an Angelesco system (see [36], 
[24], [60], [27], [28], [2]). 

In the second case, f± = /, f 2 = / 2 , where 

( 4 ) f(z) = (z 2 - l) 1/4 (z - a,y 1/2 , /(oo) = 1, a0M, 


and the pair of functions f ±, f 2 forms a (generalized) Nikishin system (see 
[43], [27], [28], [2], [5], [49], [34], [35]). 


Remark 1. In the theory of Hermite-Pade polynomials for a collection of 
three functions [1, /i, f 2 ] two opposite situations are usually selected, which 
are connected with the distribution of the branch points of the functions 
fi and f 2 . In the first case the sets of branch points A\ = A{fi) and 
A 2 = A{f 2 ) do not intersect each other. We say that the pair of functions 
forms an Angelesco system , see [28], [2], [5]. In the second case the sets of 
branch points A\ and A 2 of the functions f\ and f 2 are equivalent. We say 
that two functions fi, f 2 forms a Nikishin system , see [43], [28], [2], [5]. In 
the first case the interaction matrix M for the theory-potential-equilibrium 

/2 l\ / 2 —l\ 

vector problem is M = I ^ I. In the second case it is M = I j , see 


[25], [28], [2], [49]. The natural expectation is that such differences between 
the distribution of the branch points and the structure of the interaction 
matrices leads to essentially different limiting distributions of the zeros of 
the corresponding Hermite-Pade polynomials. However, this is not always 
the case. Figures 14, 15, 16 and 17 show the distribution of the zeros of the 
polynomials Qi 2 o,o (blue points), Qi20,i ( re d points), (Ji20,2 (black points) 
for two pair of functions: one with different branch points 


( 5 ) 


fi{z) = V( z ~ !)/(« + !). hiz) = Vi z ~ 2)/(z + 2) 
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and the other with coincident branch points 


(6) /i0) = 


z-iy / 3 

z + l 


z -2 
Z ~h 2 


1/3 


f'2 (") = 


z-iy / 3 

z + l) 


z-2 
z 2 


1/3 


From the results on figures 14-17 it follows that the distribution of the zeros 
of Hermite-Pade polynomials for these two different systems of functions will 
be equivalent (with regard to the fact that supports of the limit measures for 
the polynomials Q n ,j change over to each other). Note that the system (6) 
is a special case of a two Markov functions system, which was considered by 
E. A. Rakhmanov in [49]. The numerical results for the system (6) obtained 
here are in a good agreement with the results of E. A. Rakhmanov [49]. 


Remark 2. Instead of the pair of functions (5) with second order branch 
points we can consider a pair of functions with arbitrary order of the branch 
points. For example, the pair of functions 

(7) h(z) = y( 2 -l)/(z + l), h(z) = V(z- 2)/(z + 2) 

have the same distribution of the zeros of the corresponding Hermite-Pade 

polynomials as (5); see fig. 18-20 and comp. fig. 24-26. However, the pair 

of functions 

(8) fi(z) = \/(z- l )/(2 + l), Hz) = y(z-2)/{z + 2) 

have another distribution of the zeros of the corresponding Hermite-Pade 
polynomials, see fig. 21-23 and comp. fig. 27-29. Thus, the distribution 
depends non only of the type but also on the degree of the branch points (in 

(8) for the first function f\ instead of degree 1/3, as in (7), we took degree 

lul¬ 
lin' distribution of the zeros remains stable while moving the branch 
points of the function f\ along the imaginary axis, for example, when we 
change the points £ = ±1 into the points z = ±1 + i ■ 0.4: 

(9) fi{z) = v /{2 - fl + i-0.4))/(z - (-1 + i ■ 0.4)), 

h{z) = y( z -2)/(z + 2), 

see fig. 24-26 and comp. (7), and: 

(10) fi{z) = \f{z - (1 +* ■ 0-4))/(2 - (-1 +1-0.4)), 

k(z) = V(z-2)l(z + 2), 

see fig. 27-29 and comp. (8). This confirms that the distribution of the 
zeros of the Hermite-Pade polynomials depends not only on the geometrical 
position of the branch points, but also on the type of the branch points. 


Remark 3. If instead of (6) we consider the pair of functions 

(ID 


/i(-) 


z-l\ 1/3 fz-2\ 1/3 
z + l) \z + 2j 


h{z) 


z- 1 
z + l 


2/3 


Z-2 
z + 2 


-1/3 


then the distribution of the zeros of the corresponding Hermite-Pade poly¬ 
nomials will be different than before, see fig. 30-31. 
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1.2. A numerous conjectures about the asymptotic behavior of the zeros of 
the Hermite-Pade polynomials of both first and second kind was, in part, for¬ 
mulated in the fundamental work of J. Nuttall [46] (see also [44], [45]). These 
conjectures served as a background for some of the numerous investigations. 
Among them the prominent results of H. Stahl [55]—[59] and A. A. Gonchar- 
E. A. Rakhmanov [24], [26], first of all should be noted, and also the results 
of A. I. Aptekarev with co-authors [1], [4], [2], [5], [6], [7]. J. Nuttall’s 
exact results [44], [45], [46] about the asymptotics of the Hermite-Pade 
polynomials for a collection of m functions [/o = l,/i ,..., f m - 1], m ^ 3, 
are based mainly on the a priori assumption of the existence of an associ¬ 
ated m-sheeted Riemann surface 94 m with a canonical decomposition into 
m sheets 94^, j = 1,... , m. The sheets defined by an Abel integral of third 
kind with purely imaginary periods and logarithmic singularities of the form 
(m — 1) log z, z = z W ~ oo^ 1 ) and — log 2 , 2 = ~ oo^\ j = 2, ... ,m, and 

having also property that the first “physical” sheet 9T 1 ) is always connected 
(see [44], [46], [48]). 2 * In such case it is possible to describe the asymptotical 
behavior of the Hermite-Pade polynomials in terms connected to this RS 
(see [46], [47], [4], [2]). However, the question how to “construct” such a RS 
in the general case when m ^ 3 remains open. 2 ’ As for the present, results 
are achieved only for the case m = 3. There are two major methods to 
obtain the results of such type. The first one mainly based on the cubic 
equation (see [1], [4], [10], also [47], [41], [42]). Using this method one can 
find an explicit representation of the so-called Nuttall’s ^-function 4 (see [4], 
[10], [11]), i.e., an Abel integral of third kind, which has the level curves 
that define the needed canonical decomposition of RS 943 into three sheets . 
The other method consists of first solving the theoretical-potential extremal 
problem, connected with the existence of the so-called Nuttall condenser 
on the Riemann sphere C (see [50], [51], [38], [71]). If such a condenser 
can be found, then the RS 943 with the needed decomposition into three 
sheets is “constructed” on the base of this condenser by a specific scheme 
(see [51], [38]). Note that for now the second method can be applied only 
when the pair /i ,/2 forms a Nikishin system [51], [71]. Namely, for some 
general enough functions /i,/ 2 , which forms a complex Nikishin system, 5 
in [51] (see also [50], [38]) an existence of the Nuttall condenser, consisting 
of two non-intersecting plates and having certain symmetrical properties, is 
proved. Such a condenser is an analogue of Stahl compact, but for the case 
of Hermite-Pade polynomials when m = 3. In [38] a scheme for constructing 
a RS with three sheets and with canonical decomposition made on the base 
of an already existing Nuttall’s condenser is proposed. 

o 

In some cases, considered by Nuttall, it appears that the multi-valued analytic func¬ 
tions fi,, fm -1 continue from the first sheet onto SK m as single-valued meromor- 

phic functions and [1, /i,..., f m - 1 ] are in a sense independent. 

5 For m = 2 the existence of a corresponding hyperelliptic Riemann surface with two 
sheets follows directly from the theorems by Stahl; see [46], [62], [7], [38]. 

4 In the case of two sheeted hyperelliptic RS such function is known as Deift’s g- function; 
see [51]. 

J Under (general) Nikishin system, which contains two functions fi, fi, we understand 
such a system for which the corresponding vector theoretical-potential equilibrium problem 
is defined by a Nikishin matrix; see [2], [5], [40]. 
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As mentioned before, the numerical results are obtained here for two op¬ 
posite cases: for a pair of functions /i ,/2 creating an Angelesco system 
(see (3)) and for a pair of functions /i = /, fi = f 2 creating a (general¬ 
ized) Nikishin system (see (4)). These new numerical results give rise to 
some new conjectures about the asymptotical properties of the Hermite- 
Pade polynomials of first kind. It is well known that the Hermite-Pade 
polynomials of first and second kind are closely related, in particular, they 
are bi-orthogonal; see [46, §2, formula (2.1.9)], also the recent works [12], 
[14], [15]. It seems that the conjectures presented here can be applied also 
to the Hermite-Pade polynomials of second kind. 

For a better exposition of the new numerical results, we state at the 
beginning some well known facts about the asymptotical behavior of the 
zeros and poles of the classical Pade approximants, i.e. zeros of the Pade 
polynomials (see (2)), and also for two-point Pade approximants [20]. These 
results and their pictures (see fig. 1-11 and fig. 77-79) are needed for the 
analysis of the numerical results, connected with the behavior of the zeros 
of the Hermite-Pade polynomials. 


2. Main results 

The main empirical results obtained in the paper are presented herein. 

2.1. Angelesco system. The first part of the numerical results of the be¬ 
havior of the Hermite-Pade polynomials for a collection of three functions 
[1, fi, f 2 ] refers to the case when each of the functions f\ and /2 has a pair 
of branch points of second order and the sets of the branch points of the 
functions f\ and /2 do not intersect. Thus, the pair of functions / 1 , /2 forms 
an Angelesco system (see (3)). About the main properties of the Angelesco 
system, see [46] , also [24] . In [24] , the first results of general character about 
the convergence of the Hermite-Pade approximants of second kind for the 
collection of functions [1, / 1 ,..., / m _ 1 ], m ^ 3 were obtained, with the func¬ 
tions fi,, fm-i creating an Angelesco system and fj = dj being Markov 
functions with non-intersecting supports A j = supper,-, j = 1, 2,..., m — 1, 
along the real axis, and containing a finite number of intervals. In [24], A. 
A. Gonchar and E. A. Rakhmanov found a new property, which they called 
pushing of the support of the equilibrium measures (see also [46], [1], [2]). 
For m = 3 this method looks as follows. Let the support sets Ai,A 2 of 
the measures < 71,02 are the non-intersecting closed intervals Ai = [ai,£>i], 
A 2 = [ 02 , 62 ] of different lengths and contained in the real axis; for defi¬ 
niteness we suppose that |Ai| < |A 2 1 and that A 2 be on the right of Ai. 
It turns out that in the case when the intervals are close enough, the sup¬ 
port F 2 C A 2 of the equilibrium measure A 2 is the interval [a* 2 , 62 ], where 
a* 2 E ( 02 , 62 ) (see (12) and fig. 12). Thus, the equilibrium measures Ai,A 2 
are absolutely continuous with respect to the normalized linear Lebesgue 
measure, the density A' x (x) behaves, in the neighborhood of the points ai 
and 61 , like (x — ai ) -1 / 2 and (61 — x) -1 / 2 , respectively, and the density 
X' 2 (x) behaves, in the neighborhood of the point 02 , like (x — a ]]) 1 / 2 and like 
(62 — x ) -1 / 2 around 62 . Thus, under a specific mutual positions the smaller 
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interval pushes the support of the equilibrium measure inside the larger in¬ 
terval; this does not happen with the support of the equilibrium measure 
for a smaller interval (see fig. 12). When [ai,&i] = [—a, 0], [ 02 , 62 ] = [0,1], 
where a E (0,1), the point a% is calculated by the following formula, found 
by V. A. Kalyagin [37] (see also [46, p. 5.3, formula (5.3.18)], [1]): 


( 12 ) 


* _ (1 - a ) 3 

° 2 9(a 2 — a + 1)' 


The authors have found, by experiments, a new property called mutual 
pushing of the supports of the equilibrium measures, in the case when the 
functions f\ and fi have a pair of branch points aq, b\ E C and 02 , 62 E C, 
{ai,6i} n {02,62} = 0, located not on one line, but on two parallel lines, 
and the intervals Ai = [ 01 , 61 ] and A 2 = [ 02 , 62 ] have different lengths, i.e. 

|Ai| > |A 2 I (see fig. 42). 


Fig. 32-33. First, when these lines are far enough from each other, there 
is no collision of the support of the equilibrium measures and the supports 
of Ai and A 2 are two non-intersecting arcs, respectively 6 , see [46], [48], [2], 
The measures Ai and A 2 are absolutely continuous with respect to the length 
of the arc \dz\, and their densities A}, j = 1,2, in the neighborhoods of the 
branch points aj, bj, behave like Chebyshev measures, i.e. ~ \z — aj\~ 1 / 2 and 
~ \z— 6 j|~ 1 / 2 , respectively. This can be seen very well on figure 32, where the 
red points are the zeros of the Hermite-Pade polynomial ( 5 i 20 ,i and the black 
points are the zeros of Qi 20 , 2 - It is obvious that the extremal compact sets F\ 
and F 2 are attracted to each other, and the zeros of the polynomials Qi 20 ,j, 
which are onto Fj, are repelled from each other and from the branch points 
Uj. bj, j = 1, 2. The zeros of the polynomial Qi 2 o,o (blue points, see fig. 33) 
form a third extremal compact Fq, which separates the compact sets F\ and 
F' 2 - The distribution of the zeros of the polynomial Q n , 0 , when n —> 00 , is 
described by the third extremal measure Ao, suppAo = Fq. Thus, the fol¬ 
lowing fact is true (see [60]). Let U n (z ) := maxjlog \Q n ,\{z)\-> log \Qn, 2 {z)\}, 
z E C. Then U n is subharmonic in C. Therefore, U n (z) = —U tln (z), where 
[i n is measure, \p, n \ ^ n, U^ ln (z) = — f log | z — C| d/i n {Q) is the logarithmic 
potential with respect to (i n . For an arbitrary polynomial Q E C[^] define 
the measure 

x(Q)= £ 

C : Q(C )—0 

which counts the number of zeros of the polynomial Q. Then 
(13) lim -fi n = lim ~x(Qn,o) = A 0 

n —¥00 Ti n^-oo Tl 

(the convergence in (13) is understood as weak convergence of measures). 


6 It is well known that in the case of an Angelesco system these two arcs in a sense are 
“ attracted!’ to each other, and in the case of a Nikishin system are they “ repelled !’ from 
each other (see [2], also [46], [48], [25], [28]). 
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Fig. 34-35. Further convergence of the parallel lines along the imaginary 
axis leads to the following. While the branch points are far enough from each 
other, collision of the equilibrium measures does not occur. However, on fig. 
34 it is clearly seen, that the upper extremal compact set F\ has strongly 
curved towards the lower extremal compact set Fi. The third extremal 
compact set Fq, as before, separates Fi and F 2 . 

Fig. 36-37. Then, under further convergence of the branch points, the 
upper extremal compact set F\ has even strongly curved towards the lower 
extremal compact set i* 2 , the support of the equilibrium measure of the 
upper compact set F\ starts to break down, the second lower compact set 
almost does not change; see fig. 36. The third extremal compact set Fq, as 
before, separates the other two compact sets from each other, but now it 
touches the second compact set F 2 . 

Fig. 38—39. Finally, under certain relative positions of the pair of branch 
points, the support of the equilibrium measure of the upper extremal com¬ 
pact Fi breaks down, while the second compact set F 2 has hardly changed, 
see fig. 38. The third extremal compact set Fo, as before, separates the 
other two compact sets from each other, and touches the second compact 
set F 2 , see fig. 39. 

Fig. 40—41. Under further convergence of the pair of branch points, the 
two arcs, which are the result of the breaking of the support of the measure 
Ai, have reached the second (lower) compact set F 2 . The second compact 
set F '2 has started to change: from the total set of black points (zeros of the 
polynomial Qiso, 2 ) several points stand out, which started to form another 
component. Thus, a second component of the support of the equilibrium 
measure A 2 started to form, i.e. the support of the equilibrium measure A 2 
started breaking down on two arcs; see fig. 40. The third extremal compact 
set Fo, as before, “seeks” to separate the other two compact sets from each 
other, but now each of the compact sets F\ and F 2 has two components. It 
is clearly seen, that the compact set Fq now crosses the compact set F 2 ; see 
fig. 41. 

Fig. 42-43. Further, the two arcs, which are the result of the breaking 
of the support of the measure Ai, cross the second compact set F 2 . The 
second compact set F 2 continues to change: from the total set of black 
points (zeros of the polynomial Qiso, 2 ) even more points stand out (than 
before), which form the second component of F 2 . Thus, the forming of the 
second component of the support of the equilibrium measure A 2 continues; 
see fig. 42. The third extremal compact set Fo crosses the compact set F 2 . 
As before, it “seeks” to separate the other two compact sets Fi and F 2 from 
each other, but now each of these compact sets has by two components. It 
is clearly seen, that at the junction of the red, black and blue points appear 
two equilateral triangles with multicolored vertexes, see fig. 43. 

Fig. 44—45. The two arcs, which are the result of the breaking of the 
support of the measure Ai, even further cross the second compact set F 2 . The 
second compact set F 2 continues to change: from the total set of black points 
(zeros of the polynomial Qiso, 2 ) even more points stand out (even than 
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before), which form the second component of F 2 . Thus, the forming of the 
second component of the support of the equilibrium measure A 2 continues. 
It is clearly seen, that at the junction of the red, black and blue points 
appear two equilateral triangles with multicolored vertexes. By analogy 
with classical Pade approximants and two-point Pade approximants (see 
fig. 4 and 78), it is natural to assume, that the center of each triangle has 
a Chebotarev point v\, V 2 with zero density. At the branch points aj, bj the 
density of the measures Ai and A 2 are proportional to \z— aj|" 1//2 , \z— 6 J j~ 1 / 2 , 
j = 1,2, respectively. There is a Froissart singlet (blue) on the imaginary 
axis; see fig. 45, 44. 

Fig. 46—47. Finally, under certain relative positions of the pairs of branch 
points, the support F 2 of the equilibrium measure A 2 is separated on two 
practically equivalent arcs. However, according to the distribution of the 
zeros of the polynomial Q 180 , 2 , the density X' 2 of the equilibrium measures 
of each arc must be different. On the upper arc it behaves like a Chebyshev 
measure, that is at the end points 02,62 the density is proportional to \z — 
& 2 | —2 and |z — 62 1 -1 / 2 , respectively. The end points of the lower arc V\,V 2 
are the Chebotarev points and their density is proportional to \z — ui) 1 / 2 
and \z — U 2 I 1 / 2 ; see fig- 46. At the junction of the red, black and blue points 
appeared two equilateral triangles with multicolored vertexes, and the center 
of each has a Chebotarev point v±,V 2 ', see fig. 47. 

2.2. Nikishin system. In the theory of Pade approximants it is well known 
the property called “ Froissart doublets” , which was experimentally found 
(see [16], [30], [29]) and means that in the maximal' domain of holomorphy 
D = D(f ) of the multi-valued function / for some n G A, A C N is an infinite 
sequence, are positioned pairs of zero-pole of a diagonal Pade approximant 
(that is, different from each other zeros of the Pade polynomials P n ,o, Pn,\i 
see fig. 2-11). For every fixed n G A, these points are different from each 
other, the pole is not dependent of any singularity of the original function, 
the zero and the pole are infinitely close to each other when n —> 00 , that is 
they are asymptotically “canceled out”. In other words, the residue of the 
Pade approximant at such a pole converges to zero when n —> 00 . Because 
such poles do not correspond to the singularities of the original function, 
sometimes they are called “spurious” poles and zeros or “defects” of the 
diagonal Pade approximant. In a typical case these poles and zeros are dense 
on the Riemann surface C when n —> 00 , n G A (see [61], [64], [ 66 ], [ 68 ], 
[70]). Thus they are sometimes called wondering poles and zeros or floating 
poles and zeros. For an arbitrary algebraic function / the number of these 
pairs depends mainly on the genus of the corresponding Riemann surface, 
and also of the number of zeros of the functions A/(()) = (/ + — /“)(£), 
C G S', on the Stahl compact set S = S(/). It is shown in [ 68 ], that 
the appearance of the Froissart doublets is due to points of an “incorrect” 
interpolation of the diagonal Pade approximants [ n/n\f := —P n ,o/Pn,i in the 

~This notion has been introduced by H. Stalil [55]—[57] . The existence of such a domain 
D = D( f) = C \ S follows from the classical theorem of Stahl [62]. In this regard, the 
domain D is usually called Stahl domain and the symmetrical compact set S - Stahl 
compact for the multi-valued analytic function /. 
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Stahl domain D(f) = C \ S with another branch / of the original function 
/ when n € A. Namely, the existence of Froissart doublets does not allow 
uniform convergence of the Pade approximants in the Stahl domain (for 
details see below). 

The second part of the numerical experiments are in the case, when the 
sets of singularity points for the two functions f\ and fi intersect each 
other. Specifically, we select a collection of three functions [1,/, / 2 ], where 
the function / is of the type (4) and thus, the pair of functions /, f 2 forms 
a Nikishin system (see [5], [49]). In this case another new property has been 
found. Namely, the appearance of triple zeros ( Froissart triplets, see below), 
i.e. zeros of the Hermite-Pade polynomials Q n ,o,Qn,uQn,2, which are very 
close to each other, but still have different values. These zeros are in the 
domain of holomorphy of the functions /,/ 2 , do not correspond to either 
zeros, nor singularities of these function, for each n they are practically 
identical to each other and with the transfer from n to n + 1 they shifted in 
the complex plane as one unit. It is appropriate to compare these triplets 
with the very well known Froissart doublets for classical Pade polynomials 
(see [23], also [16], [30], [31], [29], [33], [18], [17], [62], [19]). Froissart doublets 
are sometimes called “defects” [16, Chapter 2, §2.2], and also spurious or 
wondering (floating) zeros and poles of the Pade approximant [16], [22] (see 
also [ 68 ]). 

It is considered that the existence of such zeros and poles of the Pade 
approximant does not allow uniform convergence of the Pade approximant. 
Because such zeros and poles are infinitely close to each other asymptotically, 
then when the limit is taken, they are practically “canceled out”. For some 
classes of hyperelliptic functions /, which allow the representation / = a, 

2g+2 

where the support S a = [J [e 2 j-i,& 2 j] C M of the measure a consists of 

3 = 1 

finite number of non-intersecting intervals, it was shown [65], [69], [9] that 
the movement of such poles is subject to certain regularity. Namely, the 
corresponding divisor of the Nuttall function (see [65], [69], [9]) moves 
along a Riemann surface DI 2 with two sheets of genus g ^ 1 and is subject 
to the general Dubrovin system: 


(14) 


Zk = ~ 


2w(z k ) 


n {Zk - Zj) 


L 


e 2g+2 


no *- z j 

3j£k _ 

w(x ) 


dx, 


k = l,...,g, 


where z k = dz k /dt, e. 2 g +2 is the rightmost point of the support S a = 

2g+2 

_| \e 2 j-i , e, 2 j\ C M of the measure a (it is supposed, that the support of the 
3 = 1 

measure a has g ^ 1 gaps, the endpoints of the intervals of the support are 
numbered according to the ascending values, and the path of integration in 
(14) is part of the real axis [e 2 9 + 2 ,+oo); for details about the notations see 
[66]). The points zi(i),..., z g (t), t e M+, are on the corresponding hyperel¬ 
liptic Riemann surface 912 of genus g. If r is a real-valued rational function, 
which has poles only outside S a , then for / = a + r we have: [n/n]f —> f in 
the spherical metric locally uniformly in C \S a ([n/n]f := —P n ,i/Pn,o)'i thus, 
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each pole r attracts exactly the number of poles of [n/n\f, as its multiplicity 
(see [28]). In this case the movement of the poles and the points of inter¬ 
polation of the Pade approximants [n/n\f, which are in the gaps between 
intervals, are subject to the system (14). 

The fact, that such poles of the Pade approximants form a pair with the 
zeros that are near them and in this pair they are asymptotically close to each 
other and when taking the limit are canceled out, has led some authors to 
believe that their appearance is random and is not connected with the nature 
of the original function /. Such an approach to this property reflected on 
the respective terminology: such zeros and poles of the approximant began 
to be called “spurious”. Their appearance, when using Pade approximants, 
was considered especially negatively in [16], [17], [63], therefore when n E A 
the Pade approximants [n/n]f themselves sometimes were thought to be 
defective or entirely excluded from the research [16, Chapter 2, § 2.3] or was 
proposed of using the so-called “purification” of the spurious poles [61]. 

In [32] Dumas researched the problem of the asymptotical behavior of the 
sequence {[n/n]d}neN f° r elliptic 8 functions of the special form 

(15) d(z) = \J[z — ei) • • • (z - e 4 ) - z 2 + (ei H-be 4 )z/ 2 , 

where the points ei,..., e 4 E C are pairwise different and such root branch 
is selected, that the main member of which, in a neighborhood of the point 
z = oo, is equal to z 2 ; thus d E Jf(oo). Particularly, Dumas has shown 
that in “ general position” the set of poles of the Pade approximants [n/n\d 
is dense in C. 

In [67] the result of Dumas has been extended to some classes of elliptical 
functions. It was found [68, § 1, par .2] the following property: for each 
n, from some subsequence A, always exists a point (3 = f3 n £D = C\ S, 
for which [n/n]d(/3 n ) = d(/3 n ), where d{z) = —\J{z — e 4 ) • • • (z — e 4 ) — z 2 + 
(ei + • • • + e±)z/ 2 is another branch of the elliptical function (15) in D. 
Thus, for each n, from some subsequence A, the diagonal Pade approximants 
[n/n\d interpolate at some point from the domain D another branch d of the 
function d. In the case of general position, such points of the “incorrect” 
interpolation {/ 3 n }neA are dense in C. In [68] were obtained much more 
general results in this direction. 

In the work of E. A. Rakhmanov [52] was obtained an electrostatic in¬ 
terpretation of the poles of the diagonal Pade approximant [n/n]f for some 
algebraic functions / for each fixed n. In this framework, the spurious poles 
play a special role: they should be considered a part of the “external field”. 

2.3. In the present work, experimentally was obtained a new property for 
the Nikishin system f\ = f, fi = / 2 , where 

(16) f(z) = (z 2 — l) 1 / 4 ^ — a) -1 / 2 , where /(oo) = 1 , a ^ M, 


Q 

Speaking about elliptic functions, we follow to the terminology of the monograph 
[54, Chapter 10, par. 10.10], where under elliptic function is understood a single-valued 
function defined on an elliptical Riemann surface. 
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which is connected with the behavior of the zeros of the Hermite-Pade poly¬ 
nomials of first kind. Namely, there appear, for some n E A = A(/), sin¬ 
gle spurious zeros of the Hermite-Pade polynomials, and also triple spuri¬ 
ous zeros of the Hermite-Pade polynomials. Under “spurious” zeros of the 
Hermite-Pade polynomials we understand such zeros that, first, are different 
from each other (they cannot be canceled out), second, do not correspond 
to either zeros, nor singularities of the original function, and third, signifi¬ 
cantly change their location, when we transfer from n to n + 1 (sometimes 
they may disappear). By analogue with Pade polynomials, it is naturally 
to call such zeros Froissart singlets and Froissart triplets. It is natural to 
assume the same as for Pade approximants, in the “typical” case (i.e. for 
the branch points a, which are in “general position”, that is a (j iM, see (16)) 
these spurious zeros (singlets and triplets) of the Hermite-Pade polynomials 
are dense only in “their” domain, which is the difference from Pade approx¬ 
imants. From the numerical experiments it follows that (see fig. 58, 59, 76) 
the zeros of the polynomials Q n> o and Q Ut 2 of the functions of type (16) when 
a = i ■ 1.6 have the same limiting distribution, the corresponding boundary 
compact set F^ separates the Riemann surface into three domains, two of 
which have an internal boundary arc. Thus, they have the same structure 
(see fig. 79) as in the theorem of Buslaev [20], [21] for two-point Pade ap¬ 
proximants. Remark that, the single zeros of the polynomial Q U) 2 at the 
point z ~ a corresponds to a simple pole of the function f 2 at the point 
z = a. The corresponding compact set F 2 has four Chebotarev’s points 
(three of which have zero density and one has infinite density) for the equi¬ 
librium measure, which corresponds to the limiting distribution of the zeros 
of the Hermite-Pade polynomials Q n ,o,Qn, 2 ■ It is clear, that the distribu¬ 
tion of the zeros of the Hermite-Pade polynomials Q n ,o and Q n ,2 must be 
equivalent, since with the mapping / i-a 1 // the type of the singularity of 
the original function stays the same (see (16)). 

The limiting distribution of the zeros of the polynomial Q n> 1 must be 
different. That is, it must correspond to the second compact set F\ (see fig. 
56, 57 and fig. 72-75). If now we substitute F := F 1 UF 2 , F\ n F^ / 0, then 
the complement C \F, as before, consists if three domains, but now each of 
these domains has an internal boundary arc (see fig. 79). Thus, C = FUDiU 
D 2 U H 3 . Since the pair /, / 2 forms a Nikishin system, then by the analogy 
with [50], [51], [38], [71], we can associate with the collection of functions 
[1, /, f 2 ] a Nuttall condenser. After then by the analogue with [51], we can to 
describe in terms of the condenser the limiting distribution of the zeros of the 
Hermite-Pade polynomials Q n , 0 , Qn, 1 , Qn, 2 and the corresponding Riemann 
surface with the canonical (in Nuttall’s sense) partition into three sheets. 
After that, with the help of this Riemann surface, it might be possible to 
find strong asymptotics for the Hermite-Pade polynomials. The question 
about how exactly to do this stays open; a heuristical method for solving 
this problem for functions of the type f(z) = (z — ai) Ql (z — a 2 ) a2 (z — a^) 013 , 
ai + a 2 + «3 = 0, 2aj E C \ Z, is proposed in [71]. 

Figures 64-67 show the distribution of the zeros of the Hermite-Pade 
polynomials Q n , 0 (blue points), Q n .\ (red points), Q U) 2 (black points), n = 
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121,..., 130, for the collection of three functions [1, /, f 2 ], where / is “per¬ 
turbed” with respect to the function (16): 

(17) f(z) = (z 2 - 1 ) 1/4 (z - (0.1 + iV 3 • 1.6))- 1/2 . 

It is clear, that by comparison with the unperturbed case (16), the picture 
of the distribution of the zeros is entirely the same. 

3. Concluding remarks 

Thus, the main empirical results of the paper are as follows. 

3.1. In the case when the pair of functions /i, fi forms an Angelesco system, 
where the functions /i ,/2 have the type (3), there has been found numer¬ 
ically the property of the mutual pushing of the supports of the measures, 
which are equilibrium for the extremal compact sets. 

3.2. In the case when the pair of functions /i = /, /o = f 2 forms a (general¬ 
ized) Nikishin system, where the function / has the type (4), there have been 
found numerically the properties of Froissart singlets and triplets presence. 

Thus, in the present paper there have been found numerically new prop¬ 
erties related to the behavior of the zeroes of the Hermite-Pade polynomials 
of the first kind. These numerical phenomena should be researched and 
strictly justified in future. 
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- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 1. Zeros and poles of the diagonal Pade approxi¬ 
mant [131/131]/ of the function f(z) = 1 / {{z — (—1.2 + i ■ 
0.8))(z - (0.9 + i • 1.5 ))(z - (0.5 - i • 1.2))} 1 / 3 . The zeros 
and poles of the diagonal Pade approximant [n/n\f are dis¬ 
tributed, when n —> oo, according to Stahl’s theorem [62]. 
When n = 131 the zeros and poles of the diagonal Pade 
approximant [131/131]/ of the function / are distributed ac¬ 
cordingly to the electrostatical model, by E. A. Rakhmanov 
[52], and they model the Chebotarev-Stahl compact S' 131 , 
which depends on n. Froissart doublets are not present when 
n = 131. 
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Pade approximant 



- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 2. Zeros and poles of the diagonal Pade approxi¬ 
mant [130/130]/ of the function f(z) = 1 /{{z — (—1.2 + i ■ 
0.8*))(z — (0.9 + i ■ 1.5))(z — (0.5 — i ■ 1.2))} 1 / 3 , distributed 
accordingly to the electrostatical model by E. A. Rakhmanov 
[52]. There is a Froissart doublet when n = 130 (see also fig. 
3 and fig. 4). Since the genus of the Riemann surface is 1, 
there might be at most one Froissart doublet. In full compli¬ 
ance with the Rakhmanov model [52], the Froissart doublet 
“attracts” the Stahl S-compact 5130, comp. fig. 1. 
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- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 3. The poles of the Pade approximant [130/130]/ 
approximate a Chebotarev point V130 for the 5-compact S130. 
When n —> 00 we have that v n —> v is a classical Chebotarev 
point. There is one spurious pole of the Pade approximant 
[130/130]/, it is accompanied by a spurious zero of the Pade 
approximant [130/130]/ (see fig. 4). 
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Pade approximant 
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Figure 4. The Chebotarev point should not be approxi¬ 
mated by zeros of the Pade approximant [130/130]/. There 
is one spurious zero of the Pade approximant [130/130]/, it 
is accompanied by a spurious pole of the Pade approximant 
[130/130]/ (see fig. 3). 
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Figure 5. Zeros and poles of the diagonal Pade approx- 
imants [n/n]/, n = 121,..., 130, of the function f(z) = 
1/{(~-(-1-2+u0.8))(z-(0.9+z-1.5))(,2-(0.5— z-1.2))} 1 / 3 for 
each n = 121,..., 130 are distributed in the complex plane 
accordingly to the electrostatical model by E. A. Rakhmanov 
[52]. Since the genus of the Riemann surface is 1, there 
might be at most one Froissart doublet for each n. Here 
are observed 5 Froissart doublets. In full compliance with 
the Rakhmanov model [52] the n-th Froissart doublet “at¬ 
tracts” the Stahl S'-compact S n , n E {121,..., 130}, comp, 
fig. 1. 
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Pade approximant 



- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 6. For the function f(z) = \/{{z— (—1.2+F0.8))(z— 
(0.9 + i ■ 1.5))(z — (0.5 — i ■ 1.2))} 1 / 3 poles of the Pade ap¬ 
proximates [n/n]f, n = 121,..., 130, approximate a Cheb- 
otarev point. However, the presence of Froissart doublets 
changes slightly the position of the existing points v n de¬ 
pending on n E {121,..., 130}. This is in full accordance 
with the Rakhmanov model [52], 
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- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 7. Zeros of the Pade approximants [n/n\f, n = 
121,..., 130, for the function f{z) = 1 /{{z — (—1.2 + i ■ 
0.8))(z - (0.9 + i ■ 1.5))(z - (0.5 - i ■ 1.2))} 1 / 3 . The Cheb- 
otarev point should not be approximated by zeros of the Pade 
approximant [n/n]f. 
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Pade approximant 



- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 8. Zeros and poles of the diagonal Pade approxi- 
mants [ n/n]f when n = 131 and n = 132 of the function 
f(z) = 1 /{(z - (-1.2 + i ■ 0.8)){z - (0.9 + i ■ 1.5 ))(z - (0.5 - 
i ■ 1.2))} 1 / 3 . The Froissart doublets are not present when 
n = 131 and n = 132. Therefore the S'-compacts 5i3i and 
5i32 are practically the same. 
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- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 9. Poles of the diagonal Pade approximants [ n/n\f 
when n = 131 and n = 132 of the function f(z ) = 1 /{(z — 
(-1.2 + i ■ 0.8))(z - (0.9 + i • 1.5))(z - (0.5 - i ■ 1.2))} 1 / 3 . 
The Froissart doublets are not present when n = 131 and 
n = 132. Therefore, the Chebotarev points U131 and U132, ap¬ 
proximated by the poles of the diagonal Pade approximants 
[n/n\f when n = 131 and n = 132, are very close to each 
other and are practically the same. 
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Pade approximant 



- 1.5 -1 - 0.5 0 0.5 1 1.5 2 


Figure 10. Zeros of the Pade approximants [ n/n]f when 
n = 131 and n = 132 of the function f(z ) = 

0.8))(z—(0.9+i-1.5))(z—(0.5—i-1.2))} 1 / 3 do not approximate 
the Chebotarev points U 131 and U 132 . The Froissart doublets 
are not present when n = 131 and n = 132. Therefore the 
zeros the diagonal Pade approximants [n/n\f when n = 131 
and n = 132 are practically the same. 
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-4 -2 0 2 4 

Figure 11. Zeros and poles of the diagonal Pade approxi- 
mants [103/103]/ of the function f(z ) = 1 /{(z + (4.3 + i ■ 
l.))(z — (2. + i ■ -5))(z + (2. + i • 2.))(z + (1. — i ■ 3 -)){z — (4. + 

i ■ 2 .))(z — (3. + % • 5.))} . In the limit when n oo the 

zeros and poles of the diagonal Pade approximants [n/n\f 
are distributed accordingly with Stahl’s theorem [62]. Under 
fixed n = 103 these zeros and poles are distributed in a plane, 
accordingly to the electrostatical model by Rakhmanov [52]. 
Since the genus of the Riemann surface is 1, for each n there 
might be no more than 4 Froissart doublets. Here are ob¬ 
served 3 Froissart doublets. When n = 103, in full com¬ 
pliance with the Rakhmanov model, the Froissart doublets 
“attract” the Stahl 5-compact 5 io3- 







24 


Case One, Angelesco system 


J_1_1 


1_1_1_ 

_1_1_L 


- 3 - 2-10123 

Figure 12. The distribution of the zeros of the Hermite- 
Pade polynomial Q 90,1 (red points) and ( 3 go ,2 (black points) 
for the set of three functions [l,/i,/ 2 ], where fi(z) = 
yj(z + l)/z, f 2 = y/(z- 3)/z (we select such a root branch, 
that \/T" = 1). The functions fi, f 2 are of Markov 

ty pe: fi(z) = a 1 (z), f 2 {z ) = a 2 (z), where a[(x) = 
y/(l + x)/(-x)/-K, x € (-1,0), Ai = suppcri = [-1,0], 
cr '2 (x) = 1/(3 - x)/x/i r, x € (0,3), A 2 = suppcr 2 = [0,3]. 

It is known [24], [2], that the zeros of the polynomials Q Ht i 
and Q nt2 , when taken the limit, are distributed on the sup¬ 
port of the equilibrium measures S 1 = Ai and S 2 = \ 2 and in 
accordance with their densities A] (x) and X' 2 (x). The figure 
clearly shows, that S\ = Ai, but S 2 = [^,3] C A 2 , because 
| Ai| < |A 2 |. The point a 2 G (0,3) is calculated by the for¬ 
mula (12), obtained by V. A. Kalyagin [37] (see also [46], 

[1]). The density A 7 x (x) behaves, in the neighborhoods of the 
points a± and 61 , like (x — ai ) -1 / 2 and (61 — x)^ 1 / 2 , respec¬ 
tively, and the density X 2 (x) behaves, in the neighborhood of 
the point a 2 , like (x — ) 1 / 2 and in the neighborhood of the 
point b 2 like (b 2 — x) -1//2 . 
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- 3 - 2-10123 

Figure 13. The distribution of the zeros of the Hermite- 
Pade polynomial Q 90,0 (blue points) and Qgo ,2 (black points) 
for the set of three functions [l,/i,/ 2 ], where fi(z) = 

(z + l)/z, fg = \j(z — 3) jz. The zeros of the polynomial 
<390,0 (blue points, comp fig. 33) create a third extremal 
compact Fg, which “separates” the compacts F\ and F- 2 - 
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Angelesco system and Nikishin system 
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Figure 14. The distribution of the zeros of the Hermite- 
Pade polynomials Qi20,o (blue points), Qi20,i (red points) 
and Qi20,2 (black points) for a set of three functions [1, /i, f 2 \, 
where f^z) = y/ (z - 1 )/{z + 1), f 2 (z ) = y/(z - 2 )/(z + 2). 

The sets of branch points of the functions f\ and f 2 do 
not intersect each other, and thus the pair fi,f 2 cre¬ 
ate an Angelesco system. However, the distribution of 
the zeros of the Hermite-Pade polynomials for this system 
is the same as the distribution for the Nikishin system 
AW = ((2-1)/(2+1)) 1 / 3 (( 2 -2)/(z + 2)) 1 '' 3 , f 2 (z) = 

«2 - 1 )/(z + 1)) 2/3 ((z - 2 )/(z + 2)) 1/3 ; see fig. 15-17. 











Angelesco system and Nikishin system 


27 


2 

1.5 

1 

0.5 

0 

0.5 

-1 

1.5 

-2 



-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 

Figure 15. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,o (blue points), Qi 20 ,i (red points) 
and Qi 20,2 (black points) for a set of three functions [1, /i, f 2 \, 
where h(z ) = ((z - l)/(z + 1)) 1/3 ((z - 2)/(z + 2)) 1 / 3 , 
f 2 (z) = ((z-1 )/(z + 1)) 2/3 ((z- 2)/(2; + 2)) 1/3 . The pair 
/ 11/2 create a Nikishin system, because the branch points 
of the functions are equivalent. However, the distribution of 
the zeros of the Hermite-Pade polynomials for this system is 
the same as the distribution of the zeros for the Angelesco 
system fi(z) = yj (z - 1 )/(z + 1), f 2 (z) = yj (z - 2 )/(z + 2); 
see fig. 14, 16, 17. 
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Angelesco system and Nikishin system 
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Figure 16. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,o (blue points) and Q 120,2 (black 
points) for a set of three functions [l,/i,/ 2 ], where fi(z) = 

V( z ~ l )/( z + Oi M z ) = V( z ~ 2 )/(~ + 2 )- The P air /i) h 

create an Angelesco system. However, the distribution of 
the zeros of the Hermite-Pade polynomials for this system 
is the same as the distribution for the Nikishin system 
AW = ((2-l)/(*+l)) 1/3 ((2-2)/( Z + 2)U 3 , Mz) = 

((z - l)/(z + 1)) 2,/3 ((z - 2)/(z + 2)) 123 ; comp. fig. 17. 
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Figure 17. The distribution of the zeros of the Hermite- 
Pade polynomials Q 120,0 (blue points) and Qi 20 ,i (red 
points) for a set of three functions [l,/i,/ 2 ], where 

AW = (G - i)/G + i)) 1/s ((^ - 2)/G + 2)) 1 '' 3 , h(z) = 

{(2 - l)/(2 + l)) 2 ^ ((2 - 2 )/(z + 2)) 1/,! . The pair /p /2 ere- 
ate a Nikishin system. However, the distribution of the ze¬ 
ros of the Hermite-Pade polynomials for this system is the 
same as the distribution for the Angelesco system f\(z) = 

V( z ~ l )K z + X )> M z ) = V( z ~ 2 )/(-+ 2 ); com P- 16 - 
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Figure 18. The distribution of the zeros of the Hermite- 
Pade polynomials Qi20,o (blue points), Qi20,i (red points) 
and Qi20,2 (black points) for a set of three functions [1, /i, f 2 \, 
where fi(z) = y/(z- 1 )/{z + 1), f 2 (z) = y/(z - 2 )/(z + 2). 

The sets of branch points of the functions f\ and f 2 do 
not intersect each other, and thus the pair fi,f 2 cre¬ 
ate an Angelesco system. However, the distribution of 
the zeros of the Hermite-Pade polynomials for this system 
is the same as the distribution for the Nikishin system 
AW = ((2-1)/(2+1)) 1 / 3 (( 2 -2)/(z + 2)) 1 '' 3 , f 2 ( z ) = 

((2 - 1 )/(z + 1)) 2/3 ((z - 2 )/(z + 2)) 1/3 ; see ftg. 15-17. 
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Figure 19. The distribution of the zeros of the Hermite- 
Pade polynomials Q 120,0 (blue points) for a set of three 
functions [l,/i,/ 2 ], where fi(z) = {/(z - 1 )/(z + 1), 
f 2 (z) = y/(z — 2)/ (z + 2). The pair fi, f 2 create an 
Angelesco system. However, the distribution of the ze¬ 
ros of the Hermite-Pade polynomials for this system is 
the same as the distribution for the Nikishin system 
AW = ((z-l)/(z+l)) 1/3 ((z-2)/(z + 2)) 1 '' 3 , Mz) = 

((z - 1 )/(z + 1)) 2,/3 ((z - 2 )/(z + 2)) 123 ; comp. fig. 17. 
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Figure 20. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,i (red points) for a set of three 
functions [l,/i,/ 2 ], where /i (z) = {/(z - 1 )/(z + 1), 

/ 2 (z) = y/(z — 2)/(z + 2). The pair /i,/ 2 create an 
Angelesco system. However, the distribution of the ze¬ 
ros of the Hermite-Pade polynomials for this system is 
the same as the distribution for the Nikishin system 
AW = ((z-l)/(z+l)) 1/3 ((z-2)/(z + 2)) 1 '' 3 , Mz) = 

((z - 1 )/(z + 1)) 2,/3 ((z - 2)/(z + 2)) 123 ; comp. fig. 17. 
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Figure 21. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,o (blue points), Qi 20 ,i (red points) 
and Qi 20,2 (black points) for a set of three functions [1, /i, f 2 \, 
where fi(z) = y/(z- 1 )/(z + 1), f 2 (z ) = y/(z - 2 )/(z + 2). 

The sets of branch points of the functions /i and f 2 do not 
intersect each other, and thus the pair fi, f 2 create an An¬ 
gelesco system. However, the distribution of the zeros of 
the Hermite-Pade polynomials for this system is different 
than the distribution of the pair fi(z) = y/(z — l)/(z + 1), 
f 2 (z ) = y/{z — 2 )/(z + 2). comp. fig. 18-20. This confirms 
that the distribution of the zeros of the Hermite-Pade poly¬ 
nomials depends not only on the geometrical position of the 
branch points, but also of the type of the branch points. 
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Figure 22. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,o (blue points) for a set of three func¬ 
tions [l,/i,/ 2 ], where fi{z) = y/{z- 1 )/(z + 1), f 2 (z) = 
y/(z — 2 )/(z + 2). The pair f\, f 2 create an Angelesco sys¬ 
tem. However, the distribution of the zeros of the Hermite- 
Pade polynomials for this system is different than the dis¬ 
tribution of the pair fi(z) = y/(z — 1 )/(z + 1), f 2 (z) = 
y(z-2)/(z + 2). comp. fig. 18-20. This confirms that 
the distribution of the zeros of the Hermite-Pade polynomials 
depends not only on the geometrical position of the branch 
points, but also of the type of the branch points. 
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Figure 23. The distribution of the zeros of the Hermite- 
Pade polynomials (Ji 20 ,i (red points) for a set of three func¬ 
tions [l,/i,/ 2 ], where fi{z) = y/{z- 1 )/(z + 1), f 2 (z) = 
y/(z — 2)/(z + 2). The pair f\, / 2 create an Angelesco sys¬ 
tem. However, the distribution of the zeros of the Hermite- 
Pade polynomials for this system is different than the dis¬ 
tribution of the pair fi(z) = y/(z — 1 )/(z + 1), / 2 (^) = 
y(z-2)/(z + 2). comp. fig. 18-20. This confirms that 
the distribution of the zeros of the Hermite-Pade polynomials 
depends not only on the geometrical position of the branch 
points, but also of the type of the branch points. 
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Figure 24. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,o (blue points), Qi 20 ,i (red points) 
and Qi 2 o ,2 (black points) for a set of three functions [1, /i, f 2 \, 
where fi(z) = y/(z — (1 + i ■ 0A))/(z — (—1 + i ■ 0.4)), 
f 2 (z) = y/(z — 2)/(z + 2). The set of the branch points of 
the functions f± and /2 do not intersect each other, thus the 
pair /i ,/2 create an Angelesco system; comp. 18-20. 
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Figure 25. The distribution of the zeros of the 
Hermite-Pade polynomials Qi 20 ,o (blue points) 
for a set of three functions [l,/i,/ 2 ], where 
fi{z) = y/{z~ (l + i-0.4))/(g- (-l + t-0.4)), 

fo(z) = y/(z — 2)/(z + 2). The set of the branch points of 
the functions f± and /2 do not intersect each other, thus the 
pair /i ,/2 create an Angelesco system; comp. 18-20. 
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Figure 26. The 

distribution 

of the 

zeros of 



the Hermite-Pade polynomials Qi 20 ,\ (red points) 
for a set of three functions [l,/i,/ 2 ], where 
fi{z) = $/(z- (l + i- 0-4 ))/{z - (-1 + i • 0.4)), 

f 2 (z) = y/(z — 2)/(z + 2). The set of the branch points of 
the functions f± and /2 do not intersect each other, thus the 
pair /i ,/2 create an Angelesco system; comp. 18-20. 
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Figure 27. The distribution of the zeros of the Hermite- 
Pade polynomials Qi 20 ,o (blue points), Qi 20 ,i (red points) 
and Qi 2 o ,2 (black points) for a set of three functions [1, /i, f 2 \, 
where fi{z) = y/(z — (1 + i ■ 0A))/(z — (—1 + i ■ 0.4)), 
f 2 (z) = y/(z — 2)/(z + 2). The set of the branch points of 
the functions f± and f 2 do not intersect each other, thus the 
pair /i ,/2 create an Angelesco system; comp. 21-23. 
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Figure 28. The distribution of the zeros of the 
Hermite-Pade polynomials Qi 20 ,o (blue points) 
for a set of three functions [l,/i,/ 2 ], where 
fi{z) = V(z- (l + z-0.4))/(z- (-1 +j-0.4)), 
fo(z) = (z — 2)/(z + 2). The set of the branch points of 
the functions f± and /2 do not intersect each other, thus the 
pair /i ,/2 create an Angelesco system; comp. 21-23. 
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Figure 29. The 

distribution 
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zeros of 



the Hermite-Pade polynomials Qi 20 ,\ (red points) 
for a set of three functions [l,/i,/ 2 ], where 
fi{z) = y/(z- (l + i- 0.4 ))/(z - (-1 + i • 0.4)), 

f 2 (z) = y/(z — 2)/(z + 2). The set of the branch points of 
the functions f± and /2 do not intersect each other, thus the 
pair /i ,/2 create an Angelesco system; comp. 21-23. 
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Figure 30. The distribution of the zeros of the Hermite- 
Pade polynomials Qi20,o (blue points), Qi 20 ,i (red points) 
and Qi 20,2 (black points) for a set of three functions [1, /i, f 2 \, 
where h(z) = {{z - l)/(z + 1)) 1/3 ((z - 2)/(z + 2)) 1 / 3 , 

f 2 (z) = ((z - l)/(z + 1)) 2//3 ((z - 2)/(z + 2)) -1 / 3 . Here the 
pair fit f 2 create a Nikishin system, because the branch 
points of the functions are equivalent. 
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-3-2-10123 

Figure 31. The distribution of the zeros of the Hermite- 
Pade polynomials Q 120,0 (blue points) and Qi 20 ,i (red 
points) for a set of three functions [l,/i,/ 2 ], where 

/l« = (G - 1)/G + 1)) 1/S ((■= “ 2)/G + 2)) 1 '' 3 , f 2 (z) = 

(G - 1)/G + 1)) 2/3 (G - 2)/G + 2 ))- 1 / 3 . The pair h,h 
create a Nikishin system. 
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Figure 32. The distribution of the zeros of the 
Hermite-Pade polynomial Qi 20 ,i (red points), Qi 20,2 
(black points) in the plane for two functions 

fi(z) = l/y/(z - (-1.0 + i • 0.7 ))(z - (1.0 + i • 0.7)), 
f 2 (z) = \ / . Since 

the branch points ai = —1.0 + i ■ 0.7, b\ = 1.0 + i ■ 0.7 and 
02 = —0.65 + i ■ 0 . 7,62 = 0.65 + i ■ 0.7 of the functions 
are far enough from each other, there is no collision of the 
supports of the equilibrium measures and the supports Ai 
and A 2 are two non-intersecting arcs. The measures Ai and 
A 2 are absolutely continuous with respect to the length of 
the arc \dz\, and their densities A'-, j = 1,2, behave, in the 
neighborhoods of the branch points a,j,bj, like Chebyshev 
measures, i.e. ~ \z — Oj|~ 1//2 and ~ \z — 6 J -|~ 1 / 2 , respectively. 
It is obvious that the extremal compacts F\ and F 2 are 
converging to each other, and the zeros of the polynomials 
Qi20,j, which are onto Fj, are diverging from each other and 
from the branch points a,j,bj, j = 1,2. The zeros of the 
polynomial Qi 20 ,o (blue points, see fig. 33) create a third 
extremal compact Fq, which separates the compacts F\ and 

f 2 . 
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Figure 33. The distribution of the Hermite-Pade 
polynomials Qi20,o (blue points), Qi 20 ,i (red points), 
Qi 20,2 (black points) in the plane C z for two functions 
fi(z) = l/y/{z - (-1.0 + i • 0.7))(z - (1.0 + i • 0.7)), 
fa{z) = 1/-\/(sr — (—.65 — z ■ 0.7)) [z — (.65 — i ■ 0.7)). The 
red points are distributed onto the support of the equi¬ 
librium measure Ai, when taking the limit, accordingly to 
its density A^. The black points are distributed onto the 
support of the equilibrium measure A 2 , when taking the 
limit, accordingly to its density X 2 . Same applies for the blue 
points. Both measures Ai and A 2 are absolutely continuous 
on the corresponding extremal curves F\ = supp Ai and 
F 2 = supp A 2 with respect to the length of the arc ds, 
and their densities behave like the Chebyshev measures 
\z — ctj| -1 / 2 , \z — j = 1,2 around the endpoints of 

the curves. The blue points (zeros of the polynomial Qi 20 ,o) 
separate the red points (zeros of the polynomial Qi 20 ,i) from 
the black points (zeros of the polynomial Qi 20 , 2 )- 
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Figure 34. The distribution of the Hermite-Pade polyno¬ 
mials Qi 80 ,i ( r ed points), Qiso ,2 (black points) for two func¬ 
tions fi(z ) = 1 /a J{_z — (—1.0 + i ■ 0.6)) (z — (1.0 + i ■ 0.6)), 
f- 2 (z) = 1/ y/(z — (—.65 — i ■ 0.6))(z — (.65 — i ■ 0.6)). The 
branch points have come closer to each other, however they 
are still far enough from each other and there is no collision 
of the supports of the equilibrium measures. It is clearly 
seen, that the upper extremal compact F\ has strongly bent 
towards the lower extremal compact F^. 
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Figure 35 . The distribution of the Hermite- 

Pade polynomials Qiso.o (blue points), Qiso.i (red 
points), Qi80,2 (black points) for two functions 

fi(z) = 7(1.0 - (-1.0 + i • 0.6)C)(1.0 - (1.0 + i ■ 0.6)C), 

f 2 (z) = ^(l.O - (-.65 - i ■ 0.6)C)(1.0 - (.65 - i ■ 0.6)(). 

The branch points have come closer to each other, however 
they are still far enough from each other and there is no 
collision of the supports of the equilibrium measures. It 
is clearly seen, that the upper extremal compact F\ has 
strongly bent towards the lower extremal compact F 2 . The 
third extremal compact To, as before, separates F\ and F 2 . 
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Figure 36. The distribution of the Hermite-Pade polynomi¬ 
als Qi 80 ,i (red points), Q iso ,2 (black points) for two functions 
fi(z) = l/y/(z - (-1.0 + i ■ 0.525))(z - (1.0 + i ■ 0.525)), 
f 2 (z) = 1 /yj(z - (-.65 - i ■ 0.525))(z - (.65 - i ■ 0.525)). It 
is clearly seen, that the upper extremal compact F\ has even 
strongly bent towards the lower extremal compact F 2 . The 
third extremal compact Fq. as before, separates F\ and F 2 . 
The support of the equilibrium measure of the upper ex¬ 
tremal compact F\ starts to break down, while the second 
compact F 2 has hardly changed. 
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Figure 37. The distribution of the Hermite- 
Pade polynomials Qiso.o (blue points), Qiso.i (red 
points), Qi 80,2 (black points) for two functions 
fi(z) = l/y/(z - (-1.0 + i ■ 0.525 ))(z - (1.0 + i ■ 0.525)), 
f 2 {z) = 1 /yj(z - (-.65 - i ■ 0.525 ))(z - (.65 - i ■ 0.525)). It 
is clearly seen, that the upper extremal compact F\ has 
even strongly bent towards the lower extremal compact 
F‘ 2 - The third extremal compact Fq, as before, separates 
F\ and F 2 . The support of the equilibrium measure of the 
upper extremal compact F\ starts to break down, while the 
second compact F 2 has hardly changed. The third extremal 
compact To, as before, separates the other two compacts 
from each other, but now it touches the second compact F 2 . 
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Figure 38. The distribution of the Hermite-Pade polyno¬ 
mials Qi 80 ,i ( r ed points), Qiso ,2 (black points) for two func¬ 
tions fi(z ) = \/^/(z — (—1.0 + i ■ 0.5)) (z — (1.0 + i ■ 0.5)), 
f 2 (z) = 1/ y/(z — (—.65 — z - 0.5))(z — (.65 — i ■ 0.5)). It is 
clearly seen, that under this position of the pair of branch 
points, the support of the equilibrium measure of the upper 
extremal compact F\ breaks down, while the second compact 
F ‘2 has hardly changed. The red points break down the sup¬ 
port of the equilibrium measure Ai of the compact F\ on two 


arcs. 
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Figure 39. The distribution of the Hermite- 
Pade polynomials Qiso.o (blue points), Qiso.i (red 
points), Qi 80,2 (black points) for two functions 
fi(z) = 1 /y/{z - (-1.0 + i • 0.5))(z - (1.0 + i • 0.5)), 

f 2 (z) = 1 /V( z ~ (—-65 — z - 0.5))(z — (.65 — i ■ 0.5)). It is 
clearly seen, that under this position of the pair of branch 
points, the support of the equilibrium measure of the upper 
extremal compact F\ breaks down, while the second compact 
F ‘2 has hardly changed. The third extremal compact Fq, as 
before, separates the other two compacts from each other, 
and touches the second compact F 2 
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Figure 40. The distribution of the Hermite-Pade polyno¬ 
mials Qi 80 ,i ( r ed points), Qiso ,2 (black points) for two func¬ 
tions /i (z) = 1/ y/(z — (—1.0 + i ■ 0.45))(z — (1.0 + i ■ 0.45)), 
f 2 (z) = 1/ \J{z — (—.65 — i ■ 0.45))(z — (.65 — i ■ 0.45)). It is 
clearly seen, that under this position of the pair of branch 
points, the two arcs, which are the result of the breaking 
of the support of the measure Ai, have reached the second 
(lower) compact F^. The second compact F 2 has started 
to change: from the total set of black points (zeros of the 
polynomial Qiso, 2 ) several points stand out, which started to 
form another component. Thus, a second component of the 
support of the equilibrium measure A 2 started to form, i.e. 
the support of the equilibrium measure A 2 started breaking 
down on two arcs. 
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Figure 41. The distribution of the Hermite- 
Pade polynomials Qiso.o (blue points), Qiso.i (red 
points), Qi 80,2 (black points) for two functions 
fi(z) = 1 /y/(z - (-1.0 + i • 0.45))(z - (1.0 + i ■ 0.45)), 
f 2 (z ) = 1 /yj(z — (—.65 — i ■ 0.45))(z — (.65 — % ■ 0.45)). It is 
clearly seen, that under this position of the pair of branch 
points, the two arcs, which are the result of the breaking 
of the support of the measure Ai, have reached the second 
(lower) compact F^. The second compact F 2 has started 
to change: from the total set of black points (zeros of the 
polynomial Qi 80 , 2 ) several points stand out, which started 
to form another component. Thus, a second component 
of the support of the equilibrium measure A 2 started to 
form, i.e. the support of the equilibrium measure A 2 started 
breaking down on two arcs. The third extremal compact Fq, 
as before, “seeks” to separate the other two compacts from 
each other, but now each of the compacts F\ and F 2 has by 
two components. It is clearly seen, that the compact To now 
crosses the compact TT 
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Figure 42. The distribution of the Hermite-Pade polyno¬ 
mials Qi 80 ,i ( r ed points), Qiso ,2 (black points) for two func¬ 
tions fi(z ) = \/^/(z — (—1.0 + i ■ 0.4 ))(z — (1.0 + i ■ 0.4)), 
f 2 (z) = 1/ y/(z — (—.65 — i ■ 0.4))(z — (.65 — i ■ 0.4)). It is 
clearly seen, that under this position of the pair of branch 
points, the two arcs, which are the result of the breaking of 
the support of the measure Ai, cross the second compact 
The second compact F 2 continues to change: from the total 
set of black points (zeros of the polynomial Qiso, 2 ) even more 
points stand out (than before), which form the second com¬ 
ponent of F' 2 - Thus, the forming of the second component of 
the support of the equilibrium measure A 2 continues. 








Case One, Angelesco system 


55 



1 

- 0.5 

0 

0.5 

Figure 

43. The 

distribution 

of the Hermite- 


Pade polynomials Qiso.o (blue points), Qiso.i (red 
points), Qi 80,2 (black points) for two functions 
fi(z) = 1 /y/{z - (-1.0 + i • 0A))(z - (1.0 + i • 0.4)), 

f 2 (z) = 1 /yji,z— (—.65 — i ■ 0.4))(z — (.65 — i ■ 0.4)). It is 
clearly seen, that under this position of the pair of branch 
points, the two arcs, which are the result of the breaking of 
the support of the measure Ai, cross the second compact 
F‘ 2 - The second compact F 2 continues to change: from the 
total set of black points (zeros of the polynomial Qiso, 2 ) 
even more points stand out (than before), which form the 
second component of F-^. Thus, the forming of the second 
component of the support of the equilibrium measure A 2 
continues. The third extremal compact To crosses the 
compact F 2 . As before, it “seeks” to separate the other two 
compacts F\ and F 2 from each other, but now it “fights” 
with two components. It is clearly seen, that at the junction 
of the red, black and blue points appear two “equilateral” 
triangles with multicolored vertexes. 
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Figure 44. The distribution of the Hermite-Pade polyno¬ 
mials Qi 80 ,i ( r ed points), Qiso ,2 (black points) for two func¬ 
tions fi(z ) = \/^/(z — (—1.0 + i ■ 0.3 ))(z — (1.0 + i ■ 0.3)), 
f 2 (z) = 1/ y/(z — (—.65 — i • 0.3))(z — (.65 — * • 0.3)). It is 
clearly seen, that under this position of the pair of branch 
points, the two arcs, which are the result of the breaking of 
the support of the measure Ai, even further cross the second 
compact F‘ 2 - The second compact continues to change: 
from the total set of black points (zeros of the polynomial 
Qi80,2) even more points stand out (even than before), which 
form the second component of F^- Thus, the forming of the 
second component of the support of the equilibrium measure 
A 2 continues. 
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Figure 45. The distribution of the Hermite- 
Pade polynomials Qiso.o (blue points), Qiso,i (red 
points) and Qi 80,2 (black points) for two functions 
fi(z) = 1 /y/{z - (-1.0 + i • 0.3))(z - (1.0 + i • 0.3)), 

f 2 {z) = 1 /y/{z — (—.65 — i ■ 0.3))(z — (.65 — * • 0.3)). It is 
clearly seen, that under this position of the pair of branch 
points, the two arcs, which are the result of the breaking of 
the support of the measure Ai, even further cross the second 
compact F‘ 2 - The second compact continues to change: 
from the total set of black points (zeros of the polynomial 
Qi80,2) even more points stand out (even than before), which 
form the second component of F-^. Thus, the forming of the 
second component of the support of the equilibrium measure 
A 2 continues. It is clearly seen, that at the junction of the 
red, black and blue points appear two “equilateral” triangles 
with multicolored vertexes. By analogy with classical Pade 
approximants and two-point Pade approximants (see fig. 

4 and 78) it is natural to assume, that the center of each 
triangle has a Chebotarev point v±, V 2 with zero density. 

At the branch points dj,bj the density of the measures Ai 
and A 2 are proportional to |z — a^j^ 1 / 2 , |z — 6j|~ 1//2 , j = 1, 2, 
respectively. There is a Froissart singlet (blue) on the 
imaginary axis. 
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Figure 46. The distribution of the Hermite-Pade polyno¬ 
mials Qi 80 ,i ( r ed points), Qiso ,2 (black points) for two func¬ 
tions /i (z) = 1/ y/(z — (—1.0 + i ■ 0.15))(z — (1.0 T i * 0.15)), 
f 2 (z) = 1/ \J(z — (—.65 — i ■ 0.15))(z — (.65 — i ■ 0.15)). It is 
clearly seen, that under this position of the pair of branch 
points, the support F 2 of the equilibrium measure A 2 is sep¬ 
arated on two practically equivalent arcs. However, accord¬ 
ing to the distribution of the zeros of the polynomial Qiso, 2 , 
the density A 2 of the equilibrium measures of each arc must 
be different. On the upper arc it behaves like a Chebyshev 
measure, that is at the end points 02,62 the density is pro¬ 
portional to | z — 02 1 and | z — 62 1 — , respectively. The 
end points of the lower arc v \, V 2 are the Chebotarev points 
and their density is proportional to \z — vi\ 1 ^ 2 and \z — f 2 | 1//2 . 
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Figure 47. The distribution of the Hermite- 
Pade polynomials Qiso.o (blue points), Qiso.i (red 
points), Qi 80,2 (black points) for two functions 
fi(z) = 1 iy/(z - (-1.0 + i • 0.15))(z - (1.0 + i ■ 0.15)), 
f 2 (z) = l /\/{z — (—-65 — i ■ 0.15))(z — (.65 — i ■ 0.15)). It is 
clearly seen, that under this position of the pair of branch 
points, the support F 2 of the equilibrium measure A 2 is 
separated on two practically equivalent arcs. The branch 
points 02,62 are on the upper arc, the end points of the 
lower arc v \, v 2 are the Chebotarev points and their density 
is proportional to \z — vi| 1//2 and \z — U 2 I 1 / 2 . At the junction 
of the red, black and blue points appeared two “equilateral” 
triangles with multicolored vertexes, and the center of each 
has a Chebotarev point v\. v^- 
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Figure 48. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), Q n ,i (red points), Q n , 2 
(black points), n = 61,..., 80, for a set of three functions 
[1, /, f 2 ], where f(z) = (1 — z 2 ) l ^ A ( 1 — i\J 3 • I. 62 ) -1 / 2 . The 
Riemann sphere is decomposed into 3 domains by the zeros of 
the Hermite-Pade polynomials, one of them contains the in¬ 
finity point, while the other two are symmetrical with respect 
to the imaginary axis. There is a pair of Froissart triplets in¬ 
side these two dornans for some n E {61,..., 80}, there is a 
pair of Froissart singlets in the complementary domains for 
some n E {61,..., 80}, there is one Froissart doublet on the 
negative part of the imaginary axis. This follows from the 
analysis of the next figures 49, 50, and 51. 
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Figure 49. The distribution of the zeros of the Hermite- 
Pade polynomials Q n .i (red points), n = 61,..., 80, for a set 
of three functions [1, /, / 2 ], where f(z) = (1 — z 2 ) l ^(l — iV3- 
1.6 z)” 1/2 . 
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Figure 50. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), n = 61,, 80, for a set 
of three functions [1 ,/,/ 2 ], where f(z) = (1 — z 2 ) 1 / 4 (l — iy/3- 
1.6 z)~ 1 / 2 . 
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Figure 51. The distribution of the zeros of the Hermite- 
Pade polynomials Q n , 2 (black points), n = 61,... ,80, for a 
set of three functions [1, /, / 2 ], where f(z) = (1 — z 2 ) 1 / 4 ( 1 — 
i\J 3-1.6 z)- 1 ! 2 . There is a simple zero of the polynomial Q n , 2 , 
n = 61,..., 80, on the positive part of the imaginary axis at 
the point z m a, a = i\J 3-1.6, corresponding to a simple pole 
of the function / 2 at the point 21 = ia. 
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Figure 52. The distribution of the zeros of the Hermite- 
Pade polynomials Q n , o (blue points), Q n ,i (red points), Qn ,2 
(black points), n = 121,.... 130, for a set of three functions 
[1, /, f 2 ], where f(z ) = (1 — ^ 2 ) x / 4 (l — i\J 3 • 1.6^) -1 / 2 . The 
Riemann sphere is decomposed into 3 domains by the zeros of 
the Hermite-Pade polynomials, one of them contains the in¬ 
finity point, while the other two are symmetrical with respect 
to the imaginary axis. There is a pair of Froissart triplets in¬ 
side these two domans for some n G (121,..., 130}, there is 
a pair of Froissart singlets in the complementary domains for 
some n G {121,..., 130}, there is a Froissart doublet on the 
negative part of the imaginary axis. This follows from the 
analysis of the next figures 53, 54 and 55. 
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Figure 53. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,i (red points), n = 121,..., 130, for a 
set of three functions [1, /, / 2 ], where f(z) = (1 — z 2 ) 1//4 ( 1 — 

W 3 - 1.6 z)- 1 / 2 . 
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Figure 54. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), n = 121,..., 130, for a 
set of three functions [1, /, f 2 ], where f(z) = (1 — z 2 ) 1/,4 ( 1 — 
iV3- l.Qz)- 1 / 2 . 





Case Two, Nikishin system 



- 2-1012 

Figure 55. The distribution of the zeros of the Hermite- 
Pade polynomials Q n) 2 (black points), n = 121,..., 130, for 
a set of three functions [1, /, f 2 ], where f(z) = (1 — z 2 ) 1 / 4 ^ — 
i\J 3-1.6z) -1 / 2 . There is a simple zero of the polynomial Q n , 2 , 
n = 121,..., 130, on the positive part of the imaginary axis 
at the point 2 ; « a, a = iV 3 • 1.6, corresponding to a simple 
pole of the function f 2 at the point z = ia. 
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- 3 - 2-10123 

Figure 56. The distribution of the zeros of the Hermite- 
Pade polynomials <3i65,o (blue points), Qi 65 .i (red points), 

Q 165,2 (black points) when n = 165 for a set of three func¬ 
tions [1, /, / 2 ], where f(z) = (1 — z 2 ) 1//4 ( 1 — is /3 • 1.6z) -1 / 2 . 

The Riemann sphere is decomposed into 3 domains by the 
zeros of the Hermite-Pade polynomials, one of them contains 
the infinity point, while the other two are symmetrical with 
respect to the imaginary axis. There are no Froissart zeros 
when n = 165. 
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Figure 57. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,i (red points), n = 121,..., 130, in the 
plane C 2 for a set of three functions [1,/,/ 2 ], where f(z) = 

(1 — ^ 2 ) 1 / 4 ( 1 — iyj 3 • 1 . 62 :) -1 / 2 . There are no Froissart zeros 
when n = 165. 
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Case Two, Nikishin system 



- 3 - 2-10123 

Figure 58. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), n = 121,..., 130, in the 
plane C 2 for a set of three functions [1,/,/ 2 ], where f(z ) = 

(1 — 2 2 ) 1 / 4 ( 1 — i\J 3 • 1 . 62 :) -1 / 2 . There are no Froissart zeros 
when n = 165. 
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- 3 - 2-10123 

Figure 59. The distribution of the zeros of the Hermite- 
Pade polynomials Q n , 2 (black points), n = 121,..., 130, in 
the plane C 2 for a set of three functions [1,/, / 2 ], where 
f(z) = (1 — z 2 ) 1 / 4 ( 1 — i\J 3 • 1.6z) -1 / 2 . There is a simple zero 
of the polynomial Q n , 2 , n = 121,..., 130, on the positive 
part of the imaginary axis at the point z « a, a = %\f. 3 • 1.6, 
corresponding to a simple pole of the function / 2 at the point 
z = ia. There are no Froissart zeros when n = 165. 
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Case Two, Nikishin system 



- 3 - 2-10123 

Figure 60. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), Q n ,i (red points), Q n ^ 

(black points), n = 127,128, for a set of three functions 
[1, /, f 2 ], where f(z) = (1 — z 2 ) l ^ A ( 1 — i\J 3 • 1.6^) -1 / 2 . The 
Riemann sphere is decomposed into 3 domains by the ze¬ 
ros of the Hermite-Pade polynomials, one of them contains 
the infinity point, while the other two are symmetrical with 
respect to the imaginary axis. There is a pair of Froissart 
triplets inside these two domans for some n £ {127,128}, 
there is a pair of Froissart singlets in the complementary do¬ 
mains for some n £ {127,128}, there is one Froissart doublet 
on the negative part of the imaginary axis. This follows from 
the analysis of the next figures 61, 62 and 63. 
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Figure 61. The distribution of the zeros of the Hermite- 
Pade polynomials Q n> i (red points), n = 127,128, for a set 
of three functions [1, /, / 2 ], where f(z) = (1 — z 2 ) 1 / 4 (l — iy/3- 
1.6 z)” 1/2 . 
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Case Two, Nikishin system 



- 3 - 2-10123 

Figure 62. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), n = 127,128, for a set 
of three functions [1 ,/,/ 2 ], where f(z ) = (1 — z 2 ) 1 / 4 (l — iy/3- 
1.6 z)~ 1 / 2 . 
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Figure 64. The distribution of the zeros of the Hermite- 
Pade polynomials Q n fi (blue points), Q n ,i (red points), Q n , 2 
(black points), n = 121,..., 130, for a set of three func¬ 
tions [1,/,/ 2 ], where / is a “perturbed” function: f(z) = 

(1 — z 2 ) 1 / 4 ( 1 — (0.1 + n/3-1.6)z) _1//2 . The Riemann sphere is 
decomposed into 3 domains by the zeros of the Hermite-Pade 
polynomials, one of them contains the infinity point, while 
the other two are symmetrical with respect to the imagi¬ 
nary axis. There is a pair of Froissart triplets inside these 
two domans for some n £ {121,..., 130}, there is a pair of 
Froissart singlets in the complementary domains for some 
n 6 {121,..., 130}, there is one Froissart doublet on the 
negative part of the imaginary axis. This follows from the 
analysis of the next figures 65, 66 and 67. It is clearly seen, 
that the picture of the distribution of the zeros is entirely the 
same compared to the unperturbed case. 
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Figure 65. The distribution of the zeros of the Hermite- 
Pade polynomials Q n> i (red points), n = 121,..., 130, for 
a set of three functions [1,/,/ 2 ], where / is a “perturbed” 
function: f(z) = (1 — z 2 ) 1 / 4 ( 1 — (0.1 + v v/3 • 1.6)z) -1 / 2 . 
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Figure 66. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), n = 121,..., 130, for 
a set of three functions [1,/, / 2 ], where / is a “perturbed” 
function: f(z) = (1 — z 2 ) 1//4 (l — (0.1 + i\J 3 • 1.6 )z)~ 1 / 2 . 
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Figure 67. The distribution of the zeros of the Hermite- 
Pade polynomials Q U) 2 (black points), n = 121,..., 130, for 
a set of three functions [1,/,/ 2 ], where / is a “perturbed” 
function: f{z) = (1 — 2 2 ) 1 / 4 (1 — (0.1 + iV%- 1.6)z) -1 / 2 . There 
are a simple zeros of the polynomials Q n , 2 , n = 121,, 130, 
on the positive part of the imaginary axis at the point z ~ a, 
a = 0.1 + i\J 3 • 1.6, corresponding to a simple pole of the 
function / 2 at the point 2 = ia. 
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Qn,0 ' 1 ~t~ Qn,lf “1“ Qn,2f~ 



Figure 68. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), Q n ,i (red points), Q n , 2 
(black points) when n = 166,..., 170 for a set of three func¬ 
tions [1, /, f 2 ], where f(z) = (1 — z 2 ) 1 / 4 (l — iV 3 • 0.9 z)^ 1 / 2 . 
The Riemann sphere is decomposed into 3 domains by the 
zeros of the Hermite-Pade polynomials, one of them con¬ 
tains the infinity point, while the other two are symmet¬ 
rical with respect to the imaginary axis. There is a pair 
of Froissart triplets inside these two dornans for two n of 
n E {166,..., 170}. This follows from the analysis of the 
next figures 69, 70, 71. There are no Froissart points in the 
complementary domains for n = 166,..., 170. 
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Qn,0 ' 1 ~t~ Qn,lf 1“ Qn,2f~ 



Figure 69. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,i (red points) n = 166,..., 170 for a set 
of three functions [1 ,/,/ 2 ], where f(z) = (1 — z 2 ) 1 / 4 (l — iy/3- 
0.9 z)~ 1 / 2 . 
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Figure 70. The distribution of the zeros of the Hermite- 
Pade polynomials Q n ,o (blue points), n = 166,..., 170, for a 
set of three functions [1, /, f 2 ], where f(z) = (1 — z 2 ) 1/,4 ( 1 — 
iV3- 0.9Z)- 1 / 2 . 
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Figure 71. The distribution of the zeros of the Hermite- 
Pade polynomials Q n> 2 (black points), n = 166,..., 170, for 
a set of three functions [1, /, f 2 ], where f(z) = (1 —z 2 ) 1 / 4 (l — 
i\J 3 • 0.9z)~ 1,/2 . There are simple zeros of the polynomials 
Q n , 2 , n = 166,..., 170, on the positive part of the imaginary 
axis at the point z ~ a, a = i\J 3 • 0.9, corresponding to a 
simple pole of the function f 2 at the point z = ia. 
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Figure 72. The distribution of the zeros of the Hermite- 
Pade polynomials Qi87,o (blue points), Q i 87 ,i (red points), 
Qi 87,2 (black points) when n = 187 for a set of three func¬ 
tions [1, /, f 2 ], where f(z) = (1 — z 2 ) 1 / 4 (l — iyj 3 • 0.9 z)~ 1 ! 2 . 
The Riemann sphere is decomposed into 3 domains by the 
zeros of the Hermite-Pade polynomials, one of them contains 
the infinity point, while the other two are symmetrical with 
respect to the imaginary axis. There are no Froissart points 
when n = 187. 
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Figure 73. The distribution of the zeros of the Hermite- 
Pade polynomial Q i 87 ,i (red points), when n = 187 for a set 
of three functions [1, /, / 2 ], where f(z) = (1 — z 2 ) 1 / 4 (l — iy/3- 
0.9 z)“ 1/2 . 
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Figure 74. The distribution of the zeros of the Hermite- 
Pade polynomial Qi 87 ,o (blue points), when n = 187 for a set 
of three functions [1 ,/,/ 2 ], where f(z ) = (1 — z 2 ) 1 / 4 (l — iy/3- 
0.9 z)- 1 / 2 . 
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Figure 75. The distribution of the zeros of the Hermite- 
Pade polynomial Q 137,2 (black points), when n = 187 for a 
set of three functions [1, /, / 2 ], where f(z ) = (1 — z 2 ) 1//4 (l — 
iy/ 3 • 0.9z)^ 1 / 2 . There is a simple zero of the polynomial 
Qi 87,2 on the positive part of the imaginary axis at the point 
z « a, a = iyj 3 • 0.9, corresponding to a simple pole of the 
function f 2 at the point z = ia. 






Qn, 0 


1 + Qn,lf + Qn,2f 2 



Figure 76. The distribution of the zeros of the Hermite- 
Pade polynomials Qi8o,o> Qi 80,2 (blue and black points), 
when n = 180 for a set of three functions [1,/, / 2 ], where 
f(z) = (1 — z 2 ) 1 / 4 ( 1 — iV 3 • 0.9 z)~ 1 / 2 . There is a simple zero 
of the polynomial Qiso ,2 on the positive part of the imagi¬ 
nary axis at the point z ~ o, a = W 3 • 0.9, corresponding to 
a simple pole of the function f 2 at the point z = ia. 
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Figure 77. The distribution of zeros and poles of the two- 
point Pade approximant (at points z = 0 and z = oo) of 
the function f(z) = y/(z — a)/(z — b), a = .5, b = 2. Here 
are selected two “same” branches of the function /: fo = 

^{z- .5 )/(z - 2), /oo = i J{z~ .5 )/{z - 2). The interval 
[1/2,2] is the Buslaev compact. 
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Figure 78. The distribution of zeros and poles of the two- 
point Pade approximant (at points z = 0 and z = oo) 
of the function f(z ) = ^(z — a)/(z — b ), a = .5, b = 2. 

Here are selected two “different” branches of the function 
/: fo(z) = sj(z- .5 )/(z - 2), /ooOO = -sj{z - .5 )/{z - 2). 

There is one Chebotarev point v = Vab = 1 of zero den¬ 
sity on the Buslaev compact. From this can be seen, that 
this point cannot be calculated using of the two-point Pade 
approximants, comp. 4; contrary to the case with classical 
Pade approximants (see fig. 3, 2) 
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Figure 79. The distribution of zeros and poles of the two- 
point Pade approximant of order [120/120] of the function 
f(z) = (z — ai) a (z — ( 12 )-°'; ct = 1/4, ai = 0.9 — i • 1.1 ,<22 = 
0 .1+ 1 - 0 . 2 , when two “quite different” branches are selected, 
/o = \/(z~ ai)/(z - a 2 ), foo = -\/{z-ai)/(z-a 2 ). The 
zeros (blue points) and the poles (red points) create a Buslaev 
compact. 
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